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We have performed new neutron diffuse scattering measurements in Fei- x Cr x solid solutions, 
in a concentration range 0<x<0.15, where the atomic distribution shows an inversion of the short 
range order. By optimizing the signal-background ratio, we obtain an accurate determination of the 
concentration of inversion xo =0.110(5). We determine the near neighbor atomic short range order 
parameters and pair potentials, which change sign at xq. The experimental results are compared 
with previous first principle calculations and atomistic simulations. 



INTRODUCTION 

FeCr alloys are a unique case in nature showing an 
inversion of the short range atomic order within a solid 
solution. Namely Cr atoms in a Fe matrix repel at low 
concentration, whether they attract at high concentra- 
tion This feature seems to be at the origin of the 
peculiar mechanical properties shown by FeCr alloys in 
the body centered cubic (bec) solid solution, which have 
important potential applications. There is now a consid- 
erable effort in the search of materials highly resistant 
to radiation damages, which will occur both in future 
hybrid-types reactors and in controlled fusion reactors 
[2|. FeCr alloys with body centered cubic (bec) struc- 
tures are the reference model to understand the behavior 
of ferritic FeCr steels, considered as leading candidates 
in most future nuclear energy options Q. Chromium 
substitution strongly changes the mechanical properties 
of iron, such as the swelling and formation of voids un- 
der irradiation, the radiation induced hardening, or the 
ductile- brittle transition. As a striking fact, the response 
of FeCr alloys under irradiation, as well as most of the 
mechanical properties, is highly non monotonic versus 
the Cr concentration, with a pronounced change of be- 
havior around x=0.1. 

Understanding this behavior requires a detailed knowl- 
edge of the microscopic interactions between atoms, 
which govern the occupancy of the lattice by atoms of 
different kinds. In a solid solution of a binary alloy, 
the atomic distribution is never fully random, and short 
range ordered structures are stabilized in a statistical 
way, with preferential occupations controlled by the elec- 
tronic configuration. In transition metal alloys, it has 
been known for decades that the electronic configura- 
tional energy can be written in terms of effective pair 
potentials between atoms j although the total cohe- 
sive energy cannot 0]. 

The origin of the short range order (SRO) inversion in 
FeCr comes from the electronic band structure. Intro- 
ducing Cr atoms in the Fe matrix leads to a lowering of 
the density of states at the Fermi level, due to the for- 
mation of a virtual bounded state As a result, the 
effective pair potential strongly varies with Cr concen- 



tration, and even changes sign at a given Cr content. 
The band structure calculations outline the dominant 
role played by Fe and Cr magnetism in the occurrence 
of the inversion. In FeCr, the short range order inversion 
was hrst of all predicted in the 80 les by ab initio calcula- 
tions of M. Hennion [8(, who extended to ferromagnetic 
alloys the models developed in the Coherent Potential 
Approximation (CPA) for paramagnetic transition metal 
alloys [4l6(- It was soon after observed using resistivity 
and neutron probes [1]. More recently, motivated by the 
important applications in nuclear industry, several inves- 
tigations of this effect were done both theoretically and 
experimentally. The SRO inversion was studied using 
Mossbauer spectroscopy and X-ray absorption fine 
structure (EXAFS) [10]. The mixing enthalpy was found 
to change sign around x=0.10 [HI]. Many recent calcu- 
lations of the pair potentials in FeCr were also made, 
using either atomistic simulations flil Il3j , first principle 
thermodynamical calculations l44l8l|. or ab initio per- 
turbation theories 



In this context, a precise determination of the pair 
potentials, as well as the concentration value where the 
short range order inversion occurs, is important. Neutron 
diffuse scattering is the best way to measure it, since the 
change of SRO directly affects the shape of the neutron 
cross section. This measurement is however difficult due 
to the small contrast between Fe and Cr neutron scat- 
tering lengths and to the low concentration range where 
the inversion occurs. We have performed neutron mea- 
surements in Fei- x Cr x alloys in the concentration range 
0<x<0.15. With respect to our previous measurements 
[if, the experimental set-up was optimized, resulting in 
an increase of the signal over background ratio by more 
than an order of magnitude. We could therefore measure 
samples down to very low Cr contents (x= 0.01) and scan 
the inversion region carefully. These improved measure- 
ments yield a more precise determination of the short 
range order parameters, making possible to evaluate the 
pair potentials up to the fifth neighbors, and to localize 
the critical concentration for the short range order inver- 
sion accurately. The pair potentials are compared with 
theoretical determinations. 
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EXPERIMENTAL CONDITIONS 



RESULTS AND ANALYSIS 



Fifteen polycrystalline samples were prepared with 
concentrations in the range 0<x<0.15. They were syn- 
thesized by CECM-Vitry (V), Cristaltech Grenoble (G), 
Gero-Neuhausen (L) and Cerac-Milwaukee (C), then 
shaped into cylinders of 30 mm length and 9 mm diame- 
ter. The samples homogeneity and the Cr concentrations 
were determined by chemical analysis. In the Fe-Cr sys- 
tem, the phase diagram shows a miscibility gap [22j so 
that with increasing Cr concentration, the solid solution 
starts to decompose into two bec phases, enriched in iron 
and chromium respectively. Although such decomposi- 
tion may in principle occur in the concentration range 
studied here, its kinetics is much slower than the kinetics 
of SRO, which was precisely determined by measuring the 
residual resistivity [H, [23 , 24 1 . The measurement of the 
self-diffusion coefficient, which follows an Arrhenius law 
with activation energy E=2.4 eV for a pre-exponential 
factor of 10~ 14 s, determines a suitable heat treatment. 

After an homogenization at 800°C, the samples were 
heated in a quartz tube at 520°C, a temperature where 
an equilibrium state of SRO is immediately reached. 
Then the temperature T was gradually decreased down 
to 430°C. The samples were kept at 430°C for a few 
hours and quenched into water. This procedure permits 
to reach a stable SRO state at 430°C and to preserve 
it during the quenching. For the highest concentration 
(x=0.15), where the bec solid solution may start to de- 
compose at 430° C, we studied the evolution of the SRO 
when the sample was annealed for much longer times. 

Neutron measurements were performed on the G6.1 
diffractometer of the ORPHEE reactor in Saclay, with 
a incident neutron beam of 4.73 A wavelength provided 
by a focusing graphite monochromator. The range of 
the scattering vector (O.I<K< 2.5 A -1 ), allowed measur- 
ing the diffuse scattering with a good accuracy without 
any contribution from the Bragg scattering. To decrease 
the contribution from inelastic scattering, the measure- 
ments were performed at low temperature (8 K) using a 
cryogenerator. A vertical magnetic field of 15 kOe pro- 
vided by an electromagnet was applied to saturate the 
sample in the direction perpendicular to the scattering 
plane (K_LH). Combined with zero field measurements, 
this procedure is used to separate the nuclear and mag- 
netic cross sections [25|. Both electromagnet and cryo- 
generator were placed inside a vacuum chamber, which 
decreased the environmental background by a factor 5. 
Combining this with the higher neutron flux due to high 
intensity neutron beam and focusing monochromator, the 
signal/background ratio was improved by a factor 12 with 
respect to that of ref. [lj . The neutron intensities were 
corrected for background, absorption and multiple scat- 
tering, and calibrated in absolute scale by measuring a 
vanadium standard. 



In the nuclear diffuse cross section of a FeCr dilute al- 
loy, the contribution of static lattice distortions can be 
neglected, since the atomic radii of Fe and Cr are very 
similar. The nuclear cross section then reduces to the in- 
coherent cross section and contribution from atomic short 
range order. It is shown in Fig. 1 as a function of the 
scattering vector K, for most samples in the region of the 
SRO inversion, and in Fig. 2 for samples with the highest 
Cr content x=0.15. For x =0, the nuclear cross section 
is K-independent and equal to the incoherent cross sec- 
tion of iron. For 0<x<0.10, it clearly decreases at small 
K values, showing a tendency to short range order. For 
x=0.107, the cross section is K-indcpendent again, which 
shows that the atomic distribution is random. Finally, for 
x=0.117 and above, the strong increase of the cross sec- 
tion at small K values shows the presence of short range 
clustering. 

Following the model of Cowley- Warren 26], the nu- 
clear cross section for a binary alloy Fei-zCr^ is ex- 
pressed as: 



x)(b Fe -b C rfS(K), (1) 



where crj nc is the incoherent nuclear cross section of 
the alloy: a lnc = (1-x) af^ c +x cr^' c , with af^ c = 0.427 
and (Jin C = 2.538 barn. bp e and bev are the nuclear 
coherent scattering lengths: b^ e = 0.951 and bev =0.352 
10~ 12 cm. The short range order function S(K), averaged 
for a polycrystalline alloy over all K orientations, writes: 



S(K) = 1 



>r-^ sinKR n 



n 



KR n 



(2) 



where z„, R„ and ot n are respectively the coordination 
number, radius and SRO parameter of the n th shell sur- 
rounding an atom placed at the origin. A positive (resp. 
negative) a„ value corresponds to an attractive (resp. re- 
pulsive) type of order between the central atom and the 
atoms of the same species in the shell considered. 

Due to the limited K range and to the averaged infor- 
mation given by the polycrystalline samples, the number 
of parameters must be restricted to the first neighboring 
shells. Moreover, it was not possible to determine the 
SRO parameters on close concentric shells separately. For 
a bec alloy, the shells 1 and 2 (R/a =0.866 and 1, where 
a is the lattice constant), 4 and 5 (R/a =1.66 and 1.73), 
are close to each other, and we have grouped them in the 
analysis. The averaged SRO parameter corresponding to 
the grouped (i,j) shell of coordination z^+Zj is defined as 
aij — (ziai+Zjaj)/(zi+Zj). The experimental cross sec- 
tions were fitted according to the above equations with 
the SRO parameters as fitting parameters. We obtain 
a good fit of the experimental data for all samples as 
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FIG. 1: Nuclear neutron diffuse cross section of Fei- x Cr x 
alloys. Solid lines are fits with the Cowley- Warren model, 
involving SRO parameters up to the 5 th shell. 




FIG. 2: Color on line. Nuclear neutron diffuse cross section 
of Fei-^Cr^ alloys annealed at 430° C for two annealing times: 
(a): x=0. 152(1), annealing time 18 hours; triangles: experi- 
ment; solid line: fit. (b) x=0. 151(1): annealing time 3.5 days; 
dots : experiment; dashed line: fit. Fits are performed with 
SRO parameters up to the 5 th shell. 



(32(2) mbarn.sr _1 .at _1 ), which ensures that all calibra- 
tions and corrections have been done properly. The SRO 
parameters are shown in Fig. 3 versus concentration. 
As expected the SRO parameters decrease when increas- 
ing the size of the coordination shell. The inversion of 
short range order is clearly seen on the ai2 parameter, 
which has the strongest value. At low concentrations 
an is negative and close to the curve -x/(l-x) which cor- 
responds to the maximum repulsion between Cr atoms. 
With increasing concentration, it changes sign and be- 
comes strongly positive, as expected for a tendency to 
segregation. The SRO parameter for the third atomic 
shell 0:3 is much smaller, but shows the same tenden- 
cies. In the 4-5 shells the SRO parameter is concentra- 
tion independent and close to zero, which corresponds to 
a random atomic distribution. The plot of the extrap- 
olated value of the correlation function S(0)= J2i z i a i 
versus x (Fig. 3), determines the inversion concentra- 
tion as xo=0. 110(5) for which S(0)=1. As noticed earlier 
27], S(0) must be zero in a canonical ensemble, but this 
infinitely narrow singularity in S(k) should never be ob- 
served experimentally. For the sample with the highest 
Cr content which belongs to the region of the miscibility 
gap at 430° C, we compare in Fig. 2 the nuclear cross 
section of samples annealed at long annealing times (18h 
and 3.5 days). The decomposition of the alloys in Cr- 
rich clusters starts to occur, as shown by the pronounced 
enhancement of the nuclear cross section at low K values. 
This enhancement cannot be correctly accounted for by 
the Cowley- Warren model, when only a few SRO param- 
eters are involved. It becomes more pronounced as the 
annealing time increases. 

Knowing the SRO parameters for a well-defined tem- 
perature state, one can deduce the pair interaction po- 
tentials in the Krivoglaz-Clapp-Moss mean field approach 
. In the high temperature limit, the correlation 
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function S(K) is related to the Fourier transform V(K) 
of the pair potential through the expression: 

1 



S(K) 



1 + 2x(l - x) 



V(K) 



(3) 



where To=430°C=703 K is the temperature for the 
quench-in state of the alloys. Eq. [3] can be linearized 
to obtain the pair potentials V, for the first neighbor 
shells. These pair potentials are shown in Fig. 4 versus 
concentration. As expected from the general behavior in 
transition metal alloys, for a given concentration, the pair 
potentials decrease when increasing the size of the coor- 
dination shell. Here, one notices that the three potentials 
seem to change sign at about the same concentration. 



shown by solid curves in Figs. 1 and 2. This proce- 
dure allowed an accurate determination of the SRO pa- 
rameters for all samples. For x=0, we found an abso- 
lute value of the nuclear cross section in perfect agree- 
ment with the calculated incoherent cross section of iron 



DISCUSSION 

The knowledge of the first neighbor pair potential is 
important, since it is the main ingredient to predict the 
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FIG. 3: Color on line. Variation of the SRO parameters en 
with concentration x in Fei-^Cr^ alloys, a) average parame- 
ter for the first two neighbor shells; the solid line is the curve 
for maximal repulsion between Cr nearest and next nearest 
neighbors; b) third shell; c) 4-5 shells; d) calculated value of 
the short range ordered function at k=0. 
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FIG. 5: Colour on line. Variation of the Pair potential V12 
for the first two shells with Cr concentration x; (red dots): 
this experiment; dashed line: ab initio calculation of ref. [a]; 
Inset: calculated potential deduced from ref. The solid 
line in the inset is a guide to the eye. 
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FIG. 4: Colour on line. Variation of the pair potentials Vj 
with concentration x in Fei-^Cr^ alloys. Solid lines are guides 
to the eye. 



FeCr phase diagram [22|, |33|, [34| , in a range of concen- 
tration and temperature where it cannot be measured 
at equilibrium, but which is important for nuclear ap- 
plications. The concentration range where the atomic 
distribution is almost random allows one to limit the 
segregation process induced by irradiation, keeping the 
mechanical properties induced by Cr substitution. 

In Fig. 5, we compare the variation of V12 versus con- 
centration with theoretical evaluations. In the ab initio 



calculation of [8( , the inversion for short range ordering is 
predicted for xo=0.25 instead of 0.11. In spite of this dis- 
crepancy, the quantitative agreement about the potential 
V12 is remarkable, considering that the energy of config- 
uration represents a very small part of the total cohesive 
energy. 

The Cr impurity yields a strong deformation of the 
up band, with a localized state above the Fermi level 
and a minimum of the density of states below. The oc- 
currence of a negative moment on the Cr site (namely 
antiparallel to the Fe moment) and the change of sign of 
the pair potential versus Cr content are closely related, 
both being the consequences of this deformation. Recent 
first principle calculations, using the screened general- 
ized perturbation method 19|, 2Q| confirm the dramatic 



changes of the Fermi surface topology in the majority 
spin channel of FeCr alloys, found to occur in a nar- 
row concentration range between 0.05 and 0.10. As a 
result, the first four pair potentials decrease abruptly 
with increasing Cr content. The potential Vi is found 
to change sign at x~ 0.15, and its concentration depen- 
dence mimics that of the local Cr moment. The potential 
V12, evaluated as the weighted average of the first two 
shells (Vi2=(8Vi+6V2)/14), is however larger than the 
experimental one (inset fig. 5). As argued in ref [ljj], to 
compare with experiment, one should take into account 
the influence of the magnetic state of the alloy during 
the annealing. In contrast, calculations made in the non 
magnetic or paramagnetic state 2(| 21 1 predict a much 
smoother variation of the first pair potential. 

In real space, the repulsion of diluted Cr atoms in the 
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Fe matrix may be understood as a magnetic frustration 
effect, since Cr moments are antiferromagnetic (AF) in 
pure Cr and in the bcc Cr-rich phase. In nearest neighbor 
Cr-Cr pairs, the Cr moments should be parallel to be AF 
coupled with the surrounding Fe moments, the Cr-Cr AF 
coupling would then be frustrated. Our measurements of 
the magnetic cross section of FeCr alloys (to be detailed 
later) confirm the existence of a negative moment of - 
0.8(1) fiB on the Cr site, in rather good agreement with 



previous experimental results [291 - 4311 ] and theoretical de- 
terminations [H, 15 1 . At low Cr content (x~0. 01-0.03), 
the magnetic cross section, which reflects the perturba- 
tion induced by Cr moments on neighboring Fe sites (de- 
termined from the magnetic SRO parameters) is quite 
close to that calculated for an isolated Cr moment in the 
Fe matrix 



32] 



In Fig. 6, we compare our experimental results to the 
variation of the SRO parameter a±2 calculated by eval- 
uating the mixing enthalpy. In ref.[12j, a method was 
proposed to generalize classical many body potentials, 
which could predict the inversion of the pair potential and 
change of sign in the mixing enthalpy. Further Monte- 
Carlo simulations were performed within this approach, 
and investigated the influence of the heat treatment on 
the SRO inversion, taking into account the miscibility 
gap and the formation of stable precipitates from the 
Cr-rich phase 
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At the equilibrium temperature of 
700 K, almost the equilibrium temperature of our exper- 
iment (703 K), the concentration of inversion is evaluated 
as x~ 0.12, rather close to the experimental value. An- 
other Monte-Carlo study of the thermodynamic proper- 
ties of FeCr alloys is based on a cluster expansion of the 
configurational contributions to the mixing enthalpy 13 1 . 
It yields an inversion concentration x=0.105 at T=750 K 
in very good agreement with the experimental value. We 
note however that in both cases the calculated a±2 pa- 
rameters are larger than the experimental ones, and have 
a different concentration dependency. 

To summarize, we studied the short range order inver- 
sion in Fe-Cr alloys by neutron scattering, with a much 
better accuracy than before. The concentration of inver- 
sion is found to be xo=0. 110(5). The pair potential and 
short range parameters are compared with recent calcu- 
lations. This comparison provides a stringent test of the 
calculated pair potentials, which may be useful for future 
models and predictions of the FeCr phase diagram. 

We thank S. Michaud for her participation to the ex- 
periment and data treatment, F. Ducastelle for stimu- 
lating discussions and V. Pierron-Bohnes for a critical 
reading of the manuscript. 
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